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INTRODUCTION 


Over  the  next  few  yeats,  various  kinds  of  thermal  Imaging  and  other 
electro-optical  sensing  devices  will  be  developed  and. deployed  as  inte¬ 
gral  components  of  diverse  types  of  military  and  civilian  sy8tems. 
Under  limited  visibility  conditions,  the  performance  of  such  systems  may 
be  adversely  affected.  The  point  at  which  this  degradation  becomes 
significant  depends  upon  a  number  of  factors  such  as  the  Intended 
application,  the  nature  of  low  visibility  conditions,  and  the  system  it¬ 
self.  Both  the  designer  and  the  user  should  be  aware  of  the  environ¬ 
mental  limitations  and  should  possess  the  tools  for  an  assessment  of  the 
severity  of  these  limitations.  The  purpose  of  this  report  is  to  attempt 
to  provide  such  tools. 

Dense  haze  or  fog  occurrences  in  the  lower  atmosphere  produce  sharp 
reduction  in  visibility  which  may  persist  for  hours.  The  resulting 
attenuation  of  the  electromagnetic  radiation  by  the  suspended  haze 
and/or  fog  particles  depends  upon  their  number  concentration,  refractive 
index,  and  characteristics  of  the  size  distribution.  The  extinction  (or 
conversely,  transmission)  property  of  this  polydispersion  is  highly 
wavelength  dependent  from  the  visible  through  the  infrared;  such  depen¬ 
dence  appears  to  be  strongly  related  to  the  size  range  as  well  as  the 
shape  of  the  size  distribution. 

In  the  past,  the  drop-size  data  found  in  the  literature  were  collected, 
almost  without  exception,  with  mechanical  lmpactors  and  reduced  with  the 
aid  of  a  microscope.  1  The  collection  efficiency  of  these  lmpactors 
drops  to  nearly  zero  for  particles  of  1pm  to  2pm  radius.  While  the  tiny 
particles  less  than  lim  to  2im  radius  may  play  an  insignificant  role  in 
the  mlcrophyslcs  and  dynamics  of  cloud/fog  formation,  growth,  and  dissi¬ 
pation,  they  cannot  be  lightly  neglected  in  the  optics  of  a  hazy/foggy 
environment.  Their  importance  has  been  demonstrated  by  Hindman  and 
Helmdahl.  *  As  a  result  of  missing  the  tiny  particles,  the  drop-size 
spectra  on  the  average  often  revealed  a  unimodal  distribution.  Past 
efforts  in  modeling  a  distribution  of  this  type  have  been  made  by  re¬ 
sorting  to  standard  statistical  functions’  without  regard  to  optical 
implications  of  the  distribution.  However,  with  the  advent  of  optical 
particle  counters  and  the  resulting  finer  resolution,  the  drop-size 


‘B.  J.  Mason,  1971,  The  Physics  of  Clouds,  Oxford  University  Press,  617 
pages 

*E.  E.  Hindman  II  and  0.  E.  R.  Helmdahl,  1977,  "Submicron  Haze  Droplets 
and  their  Influence  on  Visibility  in  Fog,”  preprint,  6th  Conference 
Inadvertent  and  Planned  Weather  Modification,  Am  Meteorol  Soc,  Boston, 
MA,  10-13 

*N.  H.  Fletcher,  1962,  The  Physics  of  Ra inclouds,  Cambridge  University 
Press,  386  pages 
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spectra  on  the  average  often  do  not  appear  to  be  unlmodal,  but  unmis¬ 
takably  reveal  a  second  mode  In  what  Hindman  and  He lmdah  1 1  called  the 
haze  regime. 

During  the  past  several  years,  a  significant  amount  of  drop-size  data 
has  been  collected  by  the  NATO  countries  and  the  United  States  in  an 
effort  to  establish  a  data  base  for  the  evaluation  of  various  types  of 
military  electro-optical  systems.  An  analysis  of  these  data  shows  that 
the  drop-size  distributions  on  the  average  are  usually  blmodal  and  their 
shapes  vary  considerably  with  variations  in  visibility  (e.g. ,  Low"). 
With  these  character istics  in  mind,  we  have  undertaken  to  further  ana¬ 
lyze  these  data  and  to  develop  more  adequate  haze/fog  models,  which  are 
blmodal  in  shape  and  which  reflect  the  changes  in  drop-size  spectra  with 
visibilities.  This  report  represents  the  initial  results  of  these 
investigations.  In  the  following  section,  a  review  of  the  state  of  the 
art  in  fog  modeling  together  with  a  set  of  criteria  for  the  formulation 
of  models  will  be  presented.  Next,  the  nature  of  the  data  used  and  the 
steps  in  model  development  are  discussed.  This  discussion  is  then 
followed  by  an  evaluation  of  model  outputs  and  our  concluding  remarks  in 
two  separate  sections. 


BACKGROUND 

Two  skewed  probability  density  functions  have  most  often  been  used  in 
the  past  to  model  haze/fog  drop-size  spectra,  i.e. ,  lognormal  and 
gamma.  3>  *»  1  There  are,  then,  two  variations  of  the  gamma  function.  7»* 


2E.  E.  Hindman  II  and  0.  E.  R.  Helmdahl,  1977,  "Submicron  Haze  Droplets 
and  their  Influence  on  Visibility  in  Fog,"  preprint,  6th  Conference 
Inadvertent  and  Planned  Weather  Mbdif ication.  Am  Meteotol  Soc,  Boston, 
MA,  10-13 

"R.  D.  H.  Low,  L.  D.  Duncan,  and  Y.  Y.  R.  Hsiao,  1979,  Hi ctophy steal  and 
Optical  Properties  of  California  Coastal  Fogs  at  Fort  Otd,  ASL-TR-0034, 
US  Army  Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range,  NM, 
22  pages 

*N.  H.  Fletcher,  1962,  The  Physics  of  Ralnclouds,  Cambridge  University 
Press,  386  pages 

5A.  Kh.  Khzgian,  Ed.,  1963,  Cloud  Physics,  Israel  Progran  for  Scientific 
Translations,  Jerusalem,  392  pages 

*B.  J.  Mason,  1971,  The  Physics  of  Clouds,  Oxford  University  Press,  617 
pages 

7A.  Kh.  Khrgian  and  I.  P.  Mazin,  1952,  "Distribution  of  Drops  According 
to  Size  in  Clouds,"  Trudy  Tsen  Aero  Obs,  7:56-61  (English  version) 

*D.  Deirmendjian,  1964,  "Scattering  and  Polarization  Properties  of  Water 
Clouds  and  Hazes  in  the  Visible  and  Infrared,"  Appl  Opt,  3:187-196 
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Of  Che  two,  the  fotmer  reduced  the  two  parameters  of  the  function  to 
one,  and  the  latter  added  one  parameter  to  make  it  a  three-parameter 
function.  Following  Diermend jian's  lead,  Tamplerl  and  Tomasi’  assembled 
a  large  amount  of  fog  data  from  the  literature  and  developed  statistical 
models  on  the  basis  of  these  data.  Their  results  show  that  there  can  be 
appreciable  variations  among  fog  drop-size  spectra  resulting  from  dif¬ 
ferent  formation  processes.  It  may  be  of  some  historical  interest  to 
note  that  according  to  Khrgian,  5  Levin*  first  suggested  in  1958  that  the 
size  distribution  of  droplets  obeys  the  lognormal  law  discovered  by 
Kolmogorov  for  the  size  of  gold  particles  in  placer  deposits,  but  then 
finally  settled  on  the  gamma  distribution  (1958)  on  the  basie  of  cloud- 
chamber  experiments.  Each  function  will  undoubtedly  have  its  adher¬ 
ents.  According  to  TVomey,10  there  appear  to  be  no  theoretical  grounds 
for  stating  that  haze/fog  spectra  should  necessarily  follow  any  or  all 
of  these  distribution  functions  (or  any  other  for  that  matter),  and  the 
problem  is  largely  one  of  curve  fitting. 

After  having  examined  a  large  set  of  drop-size  spectra,  Zeuv 1  1  felt  that 
most  samples  could  be  adequately  approximated  by  a  gamma  distribution. 
When  one  inspects  a  large  numb;::  of  sample  distributions  as  we  have 
done,  it  becomes  Immediately  apparent  that  no  simple  distribution  func¬ 
tion  can  reproduce  all  of  the  features  exhibited  by  measurements.  Nor 
is  it  evident  that  such  detailed  representation  is  necessary  or  desir¬ 
able,  considering  the  state  of  the  art  in  aerosol  sanpllng.  In  fact, 
having  compared  a  number  of  different  optical  particle  counters  in  a 
carefully  designed  laboratory  investigation.  Cross  and  Fenn, 1  1  showed 


*F.  Tamplerl  and  C.  Tomasi,  1976,  "Size  Distribution  Models  of  Fog  and 
Cloud  Droplets  in  Terms  of  the  Modified  Gamma  Function,"  Tellus,  28:333- 
347 

’A.  Kh.  Khrgian,  Ed.,  1963,  Cloud  Physics,  Israel  Program  for  Scientific 
Translations,  Jerusalem,  392  pages 

*L.  M.  Levin,  1958,  "Functions  to  Represent  Drop-size  Distributions  in 
Clouds,"  Izv  Geoflz  Set,  10:1211-1221  (English  version) 

1  *S.  TWomey,  1977,  Atmospheric  Aerosols,  Elsevier,  NY,  302  pages 

1  lV.  E.  Zuev,  1970,  Atmospheric  Transparency  in  the  Visible  and  in  the 
Infrared,  Israel  Pogram  for  Scientific  Translations,  Jerusalem,  2i0 
pages 

1  *T.  S.  Cross  and  R.  W.  Fenn,  Ed.,  1978,  OPAQUE  Aerosol  Counter 
Intetcompat Ison,  25  April  1977 — 4  May  1977,  AFGL-TR-7 8-0004,  USAF 
Geophysics  Laboratory,  Hanscom  AFB,  MA,  56  pages 
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that  the  optical  devices  in  cut  tent  use  fot  measuting  diop-size  distri¬ 
butions  ate  accurate  to  within  a  factor  of  2.  Since  the  obsezved  drop- 
size  data  are  typically  bimodal,  it  would  be  logical  to  fotmulate 
bimodal  distribution  functions  to  represent  them.  Of  the  two  commonly 
used  density  functions,  the  gamma  function  is  simpler  to  handle  and  its 
various  moments  can  be  readily  generated.  The  bimodal  feature  may  be 
easily  represented  by  the  weighted  sun  of  the  two  gamma  distributions. 

Since  there  is  no  theoretical  reason  for  selecting  a  priori  a  particular 
type  of  model,  a  set  of  guidelines  by  which  the  models  are  to  be  devel¬ 
oped  must  be  established.  The  guidelines  are  as  follows: 

1.  The  model  should  be  as  simple  as  practicable. 

2.  The  model  should  reproduce,  well  within  a  factor  of  2,  the  data 
and  the  derived  quantities  such  as  liquid  water  content  and  infrared 
extinction  coefficients. 

3.  The  fitting  of  the  model  to  the  measured  data  should  not  be  an 
onerous  task. 

Following  these  guidelines,  the  lognormal  distribution,  the 
one-parameter  gamma  distribution  of  Khrgian  and  Mazln,7  the  standard 
two-parameter  gamma  distribution,  and  the  three-parameter  distribution 
of  Deitmend jian *  in  two  modes  were  each  fitted  to  the  same  sets  of  data 
and  carefully  analyzed.  The  lognormal  distribution  gave  rather  incon¬ 
sistent  results,  the  one-parameter  gamma  distribution  was  consistent  but 
not  accurate  enough,  and  the  three-parameter  gamma  distribution  was 
superior  to  either  of  the  above,  but  not  enough  over  the  two-parameter 
gamma  distribution  to  warrant  the  added  burden  of  carrying  an  extra 
parameter.  Therefore,  the  bimodal  two-parameter  gamma  distribution  was 
elected. 


DATA  USED  FOR  MODEL  DEVELOPMENT 

The  data  used  for  developing  and  fitting  the  model  were  obtained  by  the 
Atmospheric  Sciences  Laboratory  (ASL)  during  an  extensive  field  experi¬ 
ment  at  Meppen,  Germany.*  Ground-based  measurements  of  haze/fog 


7A.  Kh.  Khrgian  and  I.  P.  Mazin,  1952,  "Distribution  of  Drops  According 
to  Size  in  Clouds,”  Trudy  Tsen  Aero  Obs,  7:56-61  (English  version) 

*D.  Delrmend jian,  1964,  "Scattering  and  Polarization  Properties  of  Water 
Clouds  and  Hazes  in  the  Visible  and  Infrared,"  Appl  Opt,  3:187-196 

*A  comprehensive  technical  report  on  the  Meppen  experiment  is  in 
preparation. 
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drop-size  distributions  were  made  with  a  Particle  Measuring  Systems 
FSSP-100  light  scattering  particle  measurement  device  (commonly  known  as 
Knollenberg  Counter)  which  measures  drop  sizes  from  0.25um  to  23.5pm 
radius.  The  data  from  22  February  1978  and  4  March  1978  were  reduced  to 
obtain  an  average  size  distribution  every  5  minutes  and  were  used  for 
curve  fitting.  These  two  sets  were  chosen  because  of  the  large  number 
of  drop-size  spectra  collected  (24  hours  on  22  February  and  16-1/2  hours 
on  4  March)  as  well  as  the  broad  range  of  visibilities  encountered. 
Computed  visibilities*  ranged  from  0.04  km  to  2.7  km  on  22  February  and 
from  0.02  km  to  7.5  km  on  4  Match. 

In  the  course  of  analyzing  this  huge  amount  of  data,  we  often  noted  that 
there  was  an  apparent  relationship  between  the  ranges  of  visibility  and 
the  shape  characteristics  of  a  distribution.  To  display  these  charac¬ 
teristics,  the  visibility  range  was  segmented  and  the  mean  size  distri¬ 
butions  for  the  various  segments  were  computed.  Mean  size  distributions 
for  22  February  are  shown  in  figure  1,  and  the  mean  size  distributions 
for  the  0-  to  200-meter  visibility  segment  for  several  days  are  shown  in 

figure  2.  Inspection  of  these  two  figures  not  only  reveals  the  distinct 

bimodal  character  of  the  drop-size  spectra  but  also  shows  that  the  most 

pronounced  change  with  changing  visibility  in  these  spectra  is  manifest 

in  the  second  mode. 


The  relative  importance  of  the  two  modes  can  be  investigated  by  employ¬ 
ing  some  of  the  computations  which  are  frequently  used  in  assessing  the 
degradation  of  an  electro-optical  system.  These  computations  are  liquid 
water  content  and  extinction  coefficient.  The  extinction  coefficient, 
K^,  at  a  wavelength  X  is  computed  from 


KX 


r2  7 

it  f  i  n(r )  Qext(m,r,X)  dr 
rl 


(1) 


where  n(r)  is  the  size  distribution;  Qext(m,r,X)  is  the  Mie  extinction 
efficiency  factor  at  complex  refractive  index  m,  radius  r,  and  wave¬ 
length  X;  and  r1  and  ^  are  the  range  of  radii  over  which  n(r)  is 
defined.  Liquid  water  content  W  of  the  distribution  is  given  by 


W 


yi  r^  n(r)  dr 
*2 


(2) 


♦Throughout  the  paper,  the  values  used  for  visibility  were  calculated 
from  the  drop-size  distributions. 
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with  the  implicit  assumption  that  the  density  (of  liquid  water)  is 
unity.  Now  if  we  define  K,(r)  to  be  the  partial  integral  obtained  by 
integrating  equation  (1)  from  r^  to  r  <  ^  and  let  P,(r)  =  100 
k^(r)/Kx,  then  P,(r)  is  the  percent  of  the  total  extinction  at 
wavelength  X  resulting  from  particles  of  sizes  equal  to  or  less  than 
r.  With  a  similar  definition  of  W(r)  and  Po(r),  we  have  the  percent  of 
the  total  liquid  water  content  up  to  radius  r.  Examples  of  such 
calculations  are  shown  in  figures  3  and  4  for  visibilities  of  1.95  and 
0.258  km,  respectively.  The  10  and  90  percentile  radii  from  a  number  of 
different  visibility  values  are  presented  in  table  1.  These  graphs  show 
that  the  major  contribution  to  both  extinction  and  liquid  water  comes 
from  the  second  mode  of  the  distribution. 


Fitting  the  Bimodal  Gamma  Distribution 
The  gamma  distribution  is  defined  by  the  density  function 


r(r,a,6)  - - ~7~r  e 

a!  B 

where  a  and  6  are  the  two  parameters  subject  to  a  *  1  and  B  >  0.  Pub¬ 
lished  literature  disagrees  as  to  whether  a  is  required  to  be  an  inte¬ 
ger.  Since  restricting  a  to  integer  values  simplifies  the  computations 
considerably,  we  decided  to  adopt  this  approach  in  our  model  develop¬ 
ment.  The  statistical  moments  as  given  by  the  equation 


un(a,6)  *  /  rn  r(r,a,B)  dr  ,  (4) 

o 


are 

UjCa.B)  -  B(a  +  1) 

U2(«,S)  ■  +  1)  (<*  +  2)  (5) 

lin(a,B)  -  8n(a  +  l)(a  +2)  (a  +  n)  . 
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TABLE  1.  TEN  AND  NINETY  PERCENTILE  RADII  FOR  COMPUTATION  OF 

LIQUID  WATER  CONTENT  AND  10m  EXTINCTION  COOEFFICIENT. 
DATA  IS  FROM  MEPPEN,  GERMANY,  22  FEBRUARY  1978 


10th  Percentile  90th  Percentile 

Ri  n  Ixjc  Ei  n  Tirr 


Time 

Visibility  (km) 

K10 

Lvc 

K10 

Lvc 

0005 

1.95 

0.75 

0.6 

12 

9 

0129 

0.534 

1.0 

0.7 

17.5 

18.5 

0216 

0.407 

3.5 

3.0 

15.0 

15.0 

0247 

0.242 

3.5 

2.75 

14.0 

14.5 

0334 

1.08 

1.8 

1.0 

14.0 

14.0 

0446 

0.258 

2.7 

2.2 

13.0 

12.5 

0527 

0.172 

3.0 

2.5 

12.0 

12.0 

0638 

0.287 

3.0 

2.8 

11.5 

11.5 

0751 

0.173 

3.0 

2.5 

12.5 

12.5 

0905 

0.146 

3.0 

2.2 

11.5 

11. 0 

1018 

0.529 

2.6 

2.0 

14.5 

15.0 

1229 

1.09 

2.7 

2.0 

20.8 

21.0 

1358 

0.304 

3.2 

2.9 

17.5 

18.0 

1532 

0.100 

4.5 

4.3 

18.5 

1-7.5 

1626 

0.090 

4.6 

4.5  ’ 

18.0 

19.0 

2042 

0.86 

3.9 

3.5 

13.5 

14.0 

2200 

0.112 

3.9 

3.5 

10.5 

10.5 

2301 

0.131 

4.0 

3.7 

11.6 

11.5 

r 


Figure  3.  An  example  of  cumulative  frequency  distributions  of  the 
liquid  water  content  and  the  extinction  coefficient  at  10wm 
wavelength  for  visibility  -  1.95  km.  Note  the  relatively 
high  percentage  of  sizes  below  the  lim  radius. 


Figure  4.  An  example  of  cumulative  frequency  distributions  of  the 
liquid  watet  content  and  the  extinction  coefficient  at  the 
10u&  wavelength  for  visibility  “  0.258  km.  Note  the  almost 
total  absence  of  smaller  droplets  below  lmn  radius. 
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The  values  of  a  and  8  are  readily  obtained  from  the  first  two  moments 


8  ■  **2^1  “  »*1 

and 


a 


i^/Chj  -  u*)  -  1 


(6) 

(7) 


As  mentioned  earlier,  the  data  are  typically  blmodal  and  should  be 
fitted  with  a  blmodal  function*  In  addition,  equation  (3)  represents  a 
density  function,  while  the  number  concentration  of  a  drop-size  distri¬ 
bution  is  significantly  greater  than  unity.  This  suggests  a  model  of 
the  form 


n(r)  -  Nx  +  N2  r(r,a2.B2) 


(8) 


where  Nj ,  <*j,  Bj,  N2 ,  a^,  and  B2  are  parameters  to  be  determined. 

Using  equations  (5)  and  (8),  one  obtains  a  set  of  six  simultaneous 
nonlinear  equations  of  degrees  ranging  from  2  to  12,  which  can  be  solved 
for  the  parameters  of  equation  (8).  These  equations  would  be  difficult 
to  solve  for  the  six  unknowns.  Fortunately,  a  simple  approximate  method 
for  obtaining  these  parameters  was  found.  As  will  be  seen  later,  the 
errors  resulting  from  the  approximation  must  be  Insignificant. 

An  Inspection  of  figures  1  and  2  with  the  basic  shape  of  the  gamma 
distribution  in  mind  suggests  that  the  contribution  to  n(r)  from 
comes  from  sizes  less  than  about  1.5um  radius  and  that  the 
conribution  from  IXr.c^.fiO  comes  from  radii  greater  than  1.5um.  With 
this  in  mind,  each  drop-size  distribution  was  separated  into  two  distri¬ 
butions  with  the  1.5im  radius  as  the  separation  point.  (The  1 . 5um 
separation  point  may  be  an  instrument  artifact  according  to  Pinnlck  and 
Auvermann.  1  s)  The  first  and  second  moments  and  the  number  concentra¬ 
tions  of  particles  were  computed  independently  for  each  portion.  The 
computed  number  of  particles  provides  values  for  N.  and  N2*  The  values 
of  a's  and  B's  -are  obtained  from  equations  (6)  and  (7)  in  an  iterative 
fashion.  Since  equation  (7)  will  usually  not  provide  an  integer  value 
for  a,  the  value  so  obtained  is  rounded  to  the  nearest  positive  inte¬ 
ger.  This  value  of  a  is  then  substituted  into  the  first  moment  equation 
of  (6)  to  solve  for  B. 

Frequency  distributions  of  0^  and  a2  values  ate  shown  in  table  2.  The 
table  shows  that  6  and  3  are  the  most  likely  values  for  a.  and  a2, 
respectively.  Consequently,  these  values  were  selected  as  fixed  model 


1  *R.  G.  Pinnick  and  H.  J.  Auvermann,  1979,  "Response  Characteristics  of 
Knollenbetg  Light-Scattering  Aerosol  Counters,”  J  Aerosol  Scl,  10:55-74 
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TABLE  2.  FREQUENCY  OF  OCCURRENCE  OF  VALUES  ^  AND  <*2 


Value  of  aj 

4  5  6  7  8 

22  Feb  78  1  65  98  0  0 

4  Mat  78  10  42  40  7  11 


Value  of  oj 

0  1  2  3  4  5  6 

•22  Feb  78  8  49  12  62  24  2  5 

4  Mai  78  1  1  11  22  18  14  9 


' -"^Hipa’iiii 


parameters.  The  8^  values  were  then  recomputed  by  using  i<  -  6  accord¬ 
ing  to  the  first  moment  equation  for  all  cases  where  the  original  values 
of  (*j  differed  from  6.  A  similar  recalculation  of  some  of  the  62  values 
was  also  required.  Once  these  computations  were  completed,  curve 
fitting  techniques  were  used  to  obtain  equations  relating  Np  N2,  Bp 
and  82  to  visibility. 

To  simplify  the  data  handling  procedure  and  preclude  bias  towards  the 
lower  visibility  data,  the  visibility  range  was  segmented  into  logarith¬ 
mically  equally  spaced  intervals  and  mean  values  were  computed  for  each 
Interval.  The  results  for  Np  Bp  and  $2  together  with  the  correspond¬ 
ing  equations  obtained  by  the  least  squares  curve  fitting  technique  are 
presented  in  figures  5,  6,  and  7.  These  results  indicate  that  Np  Bp 
and  $2  can  be  related  to  visibility  V  through  the  following  formulas: 


Nl  -  446/  V 


-0.106 


Bj  -  O.oi  -  0.2  log  V  , 


82*0.63  V-°*076.  (11) 

A  preliminary  equation  for  N2  was  obtained  by  the  least  squares  proce¬ 
dure,  and  simulated  drop-size  distributions  were  computed  for  a  range  of 
visibilities.  The  N2  values  were  then  modified  to  "fine  tune”  the  model 
until  visibilities  computed  from  the  simulated  distributions  agreed  with 
the  input  visibilities.  Figure  8  shows  the  results  of  such  "fine 
tuning,"  and  the  equations  for  N2  are  given  below: 


N2  -  36  V 


for  V  <  0. 5  km 


29.4  v"1*11  for  V  -  0.5  km 


Evaluation  of  Model  Outputs 


Equations  (9)  thzough  (12)  together  with  the  values  ■  6  and  ouj  -  3 
provide  the  parameters  to  be  used  in  equation  (8)  to  simulate  a  dropsize 
distribution  corresponding  to  a  given  or  an  observed  visibility.  With  a 
distribution  known,  it  is  a  simple  matter  to  compute  the  resulting 
extinction  or  transmission  of  any  desired  wavelength.  As  discussed 
earlier,  the  model  is  not  intended  to  reproduce  all  the  details  of  the 
measurements;  the  goal  was  to  simulate  their  principal  characteristics 
and  extinction  properties.  An  example  of  the  model's  capability  to 
satisfy  the  first  objective  is  shown  in  figure  9.  The  comparison  of  the 
measurement  with  simulation  can  be  considered  excellent  for  radii 
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VISIBILITY  (km) 


5.  The  number  of  drops  la  the  fltst  mode  as  a  function  of  vial 
billty  in  a  least-squares  fit. 


6.  The  gaaaa  distribution  parameter  of  the  first  mode  as 
function  of  visibility  in  a  least-squares  fit. 


Figure  7.  the  gamma  distribution  parameter  of  the  second  mode  as  a 
function  of  visibility  in  a  least-squares  fit. 


Figure  8.  the  number  of  drops  in  the  second  mode  as  a  function  of 
visibility  in  a  fine-tuned  least  squares  fit. 
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Figure  9.  An  example  of  how  well  the  model  distribution  matches  the 
observed  dlstzlbutlon  fot  visibility  “  0.422  km. 
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smaller  than  about  10pm.  The  model  cannot  reproduce  the  apparent  third 
mode  between  10pm  and  15pm  radii.  This  apparent  third  mode,  which  can 
also  be  seen  in  figure  1,  appears  in  about  one-half  of  the  observa¬ 
tions.  This  occurrence  is  not  at  all  surprising,  as  may  be  deduced  from 
a  theoretical  study  by  Neiburger  and  Chien. 1  *  In  this  study,  they 
succeeded  in  generating  a  bimodal  drop-size  distribution  in  the  fog 
sector  after  some  50  minutes  of  droplet  growth  under  the  influence  of 
droplet  collision  and  coalescence. 

Liquid  water  content  and  extinction  coefficient  for  wavelengths  of  3.8pm 
and  10pm  computed  from  the  model  drop-size  distributions  were  compared 
with  those  derived  from  corresponding  measured  ones  of  22  February  1978 
and  4  March  1978.  As  before,  mean  values  for  specific  visibility  inter¬ 
vals  were  also  computed.  These  values  are  tabulated  in  table  3.  A 
cursory  comparison  of  the  values  of  the  corresponding  parameters  shows 
that  deviations  of  the  model  values  are  on  the  average  no  greater  than 
±15%  of  those  derived  from  observations — well  within  the  measurement 
capability  of  the  particle  counter  in  current  use. 

Concluding  Remarks  and  Some  Observations 

The  methodology  which  has  been  developed  and  exploited  in  this  report 
provides  a  capability  to  generate  bimodal  drop-size  distributions  con¬ 
sisting  of  both  the  haze  and  fog  sectors  commensurate  with  measurement 
data  (at  least  in  the  case  of  the  Meppen  fogs)  and  possesses  the  addi¬ 
tional  capability  to  provide  a  distribution  corresponding  to  a  given  or 
observed  visibility.  The  model  distribution  so  established  may  then  be 
used  to  calculate  extinction  coefficients  at  other  wavelengths. 

Although  the  equations  presented  herein  do  not  provide  the  parameters 
for  depicting  all  types  of  haze/fog  drop-size  distributions,  they  are 
nevertheless  quite  adequate  in  representing  most  of  the  observations 
taken  at  Meppen,  with  the  possible  exception  of  those  of  19  February 
1978.  It  may  thus  be  said  that  the  approach  used  for  modeling  the 
Meppen  fogs  with  a  bimodal  gamma  distribution  could  be  readily  applied 
to  other  fog  data.  In  concluding  this  report,  we  make  a  few  observa¬ 
tions  about  the  difficulties  encountered  in  haze/fog  modeling  as  a 
result  of  the  • experiences  gained  in  our  thoroughgoing  analysis  of  the 
Meppen  fogs. 

1.  There  are  strong  indications  that,  because  of  the  continual 
presence  of  submicron-sized  pollution  particles  which  may  be  too  tiny  or 
may  not  be  hygroscopic  enough  to  grow  into  fog  droplets,  there  will 
always  be  a  haze  regime  in  a  fog,  the  magnitude  of  which  depends  upon 


1  %M.  Neiburger  and  C.  W.  Chien,  1960,  "Computations  of  the  Growth  of 
Cloud  Drops  by  Condensation  Using  an  Electronic  Digital  Computer,”  Phy s 
of  Precipitation,  Geophys.  Monor.  No.  5,  Amer.  Geophys.  Union, 
Washington,  DC,  191-210 
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TABLE  3.  COMPARISON  OF  DERIVED  QUANTITIES  FROM 

MEASURED  AND  MODELED  SIZE  DISTRIBUTIONS. 
PRIMED  QUANTITIES  INDICATE  MODELED  VALUES. 


Visibility 

Lwc 

Lwc ' 

k3.8 

k3.8 

K10 

K10 

4  Match  1978 

6.640 

0.0004 

0.0003 

0.162 

0.241 

0.028 

0.026 

1.000 

0.0046 

0.0045 

3.58 

3.27 

0.512 

0.655 

0.614 

0.0087 

0.0102 

6.86 

6.43 

1.028 

1.327 

0.380 

0.0158 

0.0143 

12.79 

11.83 

1.98 

2.54 

0.201 

0.0410 

0.0440 

28.45 

26.47 

5.64 

5.02 

0.095 

0.1110 

0.0985 

56.57 

56.82 

16.74 

14.02 

0.065 

0.1750 

0.1567 

81.68 

87.81 

27.48 

22.62 

22_ 

Febiuaty 

1978 

2.300 

0.0016 

0.0018 

0.633 

1.01 

0.189 

0.158 

1.490 

0.0033 

0.0031 

1.395 

1.85 

0.453 

0.364 

0.980 

0.0065 

0.0054 

2.67 

3.67 

0.842 

0.979 

0.630 

0.0131 

0.0101 

5.27 

6.37 

2.02 

1.31 

0.390 

0.0287 

0.0180 

9.71 

10.05 

4.83 

2.42 

0.251 

0.0469 

0.0313 

17.01 

19.18 

7.64 

4.29 

0.160 

0.0691 

0.0533 

29.14 

31.87 

11.03 

7.45 

0.099 

0.1047 

0.0949 

52.59 

54.88 

16.29 

13.50 

0.064 

0.2335 

0.1599 

72.33 

89.42 

38.81 

23.09 
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the  state  of  pollution  at  a  locality.  Therefore,  one  would  expect  to 
find  bimodal  drop-size  spectra  in  a  fog  observed  with  an  optical  parti¬ 
cle  counter  capable  of  measuring  submict on-sized  particles.  As  the  fog 
evolves  and  grows,  a  third  mode  generally  appears  under  the  influence  of 
droplet  collision  and  coalescence.  Further  evidence  of  this  may  be 
found  in  Jiusto.  15  Since  droplet  collision,  collection,  and  coalescence 
are  stochastic  in  nature,  16  fog  drop-size  spectra  may  well  be  multimodal 
as  noted  by  Eldtidge.  17  Our  present  approach  merely  strikes  for  a 
compromise. 

2.  In  the  absence  of  condensation  nuclei,  there  would  be  no  haze  or 
fog  at  the  kind  of  saturation  and  supersatur ation  found  in  the  atmo¬ 
sphere.  An  abundance  of  nuclei  of  substantial  sizes  would  form  fog  in  a 
rather  slightly  supersaturated  or  barely  saturated  atmosphere,  but  the 
nuclei  are  subject  to  gravitational  settling  as  well  as  washout.  Thus, 
as  a  fog  persists  under  a  given  synoptic  situation  for  2  or  more  days, 
one  would  expect  that  its  drop-size  spectra  would  change  from  day  to  day 
for  the  above  reasons.  This  tendency  to  change  may  explain  why  the 
spectra  of  19  February  1978,  as  shown  in  figure  2,  differed  considerably 
from  the  others. 

3.  The  above  observations  then  lead  to  a  third  one.  Regardless  of 
whether  meteorological  factors  will  or  will  not  enter  into  our  modeling 
considerations,  it  is  Imperative  to  have  some  knowledge  of  the  pollution 
state  such  as  light,  moderate,  or  heavy  pollution  as  well  as  of  the 
abundance  of  cloud /fog  condensation  nuclei  at  a  locality.  With  this 
knowledge,  it  would  then  become  possible  to  refine  our  model  categories 
to  be  suitable  for  use  in  a  realistic  battlefield  environment. 


1 5J.  E.  Jiusto,  1979,  Considerations  in  the  Optical  Characterization  of 
the  Atmosphere,  ASL-CR-7 9-0100-3,  US  Army  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  NM,  26  pages 

1  ^1.  R.  Pruppacher  and  J.  D.  Klett,  1978,  Mlcrophysics  of  Clouds  and 
Precipitation,  Reldel,  Boston,  MA,  714  pages 

1 7R.  G.  Eldtidge,  1966,  "Haze  and  Fog  Aerosol  Distributions,”  J  Atmos 
Sci,  23:605-613 
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US  Arnty  Intelligence  Center  4  School 
ATTN:  ATSI-CS-C  {2LT  Coffman) 

Fort  Huachuca,  AZ  85613 

Commander 

US  Army  Yuma  Proving  Ground 
ATTN:  STEYP-MSA-TL 
Bldg  2105 
Yuma,  AZ  85364 

Northrop  Corporation 
Electro-Mechanical  Division 
ATTN:  Dr.  Richard  D.  Tooley 
500  East  Orangethorpe  Avenue 
Anaheim,  CA  92801 

Commander 

Naval  Weapons  Center 

ATTN:  Code  3918  (Dr.  Alexis  Shlanta) 

China  Lake,  CA  93555 

Hughes  Helicopters 

Army  Advanced  Attack  Helicopter  Weapons 
ATTN:  Mr.  Charles  R.  Hill 
Centinela  and  Teale  Streets 
Bldg  305,  MS  T-73A 
Culter  City,  CA  90230 

Commander 

US  Army  Combat  Developments 
Experimentation  Command 
ATTN:  ATEC-PL-M  (Mr.  Gary  G.  Love) 

Fort  Ord,  CA  93941 

SRI  International 
ATTN:  K2060/Dr.,  Edward  E.  Uthe 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

SRI  International 
ATTN:  Mr.  J.  E.  Van  der  Laan 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

Joane  May 

Naval  Environmental  Prediction 
Research  Facility  (NEPRF) 

ATTN:  Library 
Monterey,  CA  93940 


Syl vanla  Systems  Group, 

Western  Division 
GTE  Products  Corporation 
ATTN:  Technical  Reports  Library 
P.0.  Box  205 

Mountain  View,  CA  94042 

Syl vania  Systems  Group 
Western  Division 
GTE  Products  Corporation 
ATTN:  Mr.  Lee  W.  Carrier 
P.0.  Box  188 

Mountain  View,  CA  94042 

Pacific  Missile  Test  Center 
Geophysics  Division 
ATTN:  Code  3253 
Point  Mugu,  CA  93042 

Pacific  Missile  Test  Center 
Geophysics  Division 
ATTN:  Code  3253  (Terry  E.  Battallno) 
Point  Mugu,  CA  93042 

Effects  Technology  Inc. 

ATTN:  Mr.  John  D.  Carlyle 
5383  Hollister  Avenue 
Santa  Barbara,  CA  93111 

Commander 

Naval  Ocean  Systems  Center 

ATTN:  Code  532  (Dr.  Juergen  Richter) 

San  Diego,  CA  92152 

Commander 

Naval  Ocean  Systems  Center 

ATTN:  Code  5322  (Mr.  Herbert  G.  Hughes) 

San  Diego,  CA  92152 

Commander 

Naval  Ocean  Systems  Center 
ATTN:  Code  4473  (Tech  Library) 

San  Diego,  CA  92152 

The  RAND  Corporation 
ATTN:  Ralph  Huschke 
1700  Main  Street 
Santa  Monica,  CA  90406 
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Particle  Measurinq  Systems,  Inc. 

ATTN:  Dr.  Robert  G.  Knollenberg 
1855  South  57th  Court 
Boulder,  CO  80301 

US  Department  of  Commerce 
National  Oceanic  and  Atmospheric  Admin 
Environmental  Research  Laboratories 
ATTN:  Library,  R-51,  Technical  Reports 
325  Broadway 
Boulder,  CO  80303 

US  Department  of  Commerce 
National  Oceanic  and  Atmospheric  Admin 
Environmental  Research  Laboratories 
ATTN:  R45X3  (Dr.  Vernon  E.  Derr) 
Boulder,  CO  80303 

US  Department  of  Commerce 
National  Telecommunications  and 
Information  Administration 
Institute  for  Telecommunication  Sciences 
ATTN:  Code  1-3426  (Dr.  Hans  J.  Liebe) 
Boulder,  CO  80303 

AFATL/DLODL 
Technical  Library 
Eglin  AFB,  FL  32542 

Commanding  Officer 
Naval  Training  Equipment  Center 
ATTN:  Technical  Information  Center 
Orlando,  FL  32813 

Georgia  Institute  of  Technology 
Engineering  Experiment  Station 
ATTN:  Dr.  Robert  W.  McMillan 
Atlanta,  GA  30332 

Georgia  Institute  of  Technology 
Engineering  Experiment  Station 
ATTN:  Dr.  James  C.  Wiltse 
Atlanta,  GA  30332 

Commandant 

US  Army  Infantry  Center 

ATTN:  ATSH-CD-MS-E  (Mr.  Robert  McKenna) 

Fort  Benning,  GA  31805 


Commander 

US  Army  Signal  Center  &  Fort  Gordon 

ATTN:  ATZHCD-CS 

Fort  Gordon,  GA  30905 

Commander 

US  Arny  Signal  Center  &  Fort  Gordon 

ATTN:  ATZHCD-0 

Fort  Gordon,  GA  30905 

USAFETAC/DNE 

ATTN:  Mr.  Charles  Glauber 
Scott  AFB,  IL  62225 

Commander 

Air  Weather  Service 

ATTN:  AWS/DNDP  (LTC  Kit  G.  Cottrell) 

Scott  AFB,  IL  62225 

Commander 

Air  Weather  Service 

ATTN:  AWS/DOOF  (MAJ  Robert  Wright) 

Scott  AFB,  IL  62225 

Commander 

US  Army  Combined  Arms  Center 

5  Ft.  Leavenworth 

ATTN:  ATZLCA-CAA-Q  (Mr.  H.  Kent  Pickett) 
Fort  Leavenworth,  KS  66027 

Commander 

US  Army  Combined  Arms  Center 

6  Ft.  Leavenworth 

ATTN:  ATZLCA-SAN  (Robert  DeKinder,  Jr.) 
Fort  Leavenworth,  KS  66027 

Commander 

US  Army  Combined  Arms  Center 
&  Ft.  Leavenworth 

ATTN:  ATZLCA-SAN  (Mr.  Kent  I.  Johnson) 
Fort  Leavenworth,  KS  66027 

Commander 

US  Army  Combined  Arms  Center 
4  Ft.  Leavenworth 

ATTN:  ATZLCA-WE  (LTC  Darrell  Holland) 
Fort  Leavenworth,  KS  66027 
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President 

USAARENBD 

ATTN:  ATZK-AE-TA  (Dr.  Charles  R.  Leake) 
Fort  Knox,  KY  40121 

Commander 

US  Army  Armor  Center  and  Fort  Knox 
ATTN:  ATZK-CD-MS 
Fort  Knox,  KY  40121 

Commander 

US  Army  Armor  Center  and  Fort  Knox 
ATTN:  ATZK-CD-SD 
Fort  Knox,  KY  40121 

Aerodyne  Research  Inc. 

ATTN:  Dr.  John  F.  Ebersole 
Crosby  Drive 
Bedford,  MA  01730 

Commander 

Air  Force  Geophysics  Laboratory 
ATTN:  OPA  (Dr.  Robert  W.  Fenn) 

Hanscom  AFB,  MA  01731 

Commander 

Air  Force  Geophysics  Laboratory 
ATTN:  OPI  (Dr.  Robert  A.  McClatchey) 
Hanscom  AFB,  MA  01731 

Massachusetts  Institute  of  Technology 

Lincoln  Laboratory 

ATTN:  Dr.  T.  J.  Gobllck,  B-370 

P.0.  Box  73 

Lexington,  MA  02173 

Massachusetts  Institute  of  Technology 

Lincoln  Laboratory 

ATTN:  Dr.  Michael  Gruber 

P.0.  Box  73 

Lexington,  MA  02173 

Raytheon  Company 

Equipment  Division 

ATTN:  Dr.  Charles  M.  Sonnenscheln 

430  Boston  Post  Road 

Wayland,  MA  01778 


Commander 

US  Amy  Ballistic  Research  Laboratory/ 
ARRADCOM 

ATTN:  DRDAR-BLB  (Mr.  Richard  McGee) 
Aberdeen  Proving  Ground,  MD  21005 

Commander/Di  rector 
Chemical  Systems  Laboratory 
US  Army  Armament  Research 
&  Development  Command 
ATTN:  DRDAR-CLB-PS  (Dr.  Edward  Stuebing) 
Aberdeen  Proving  Ground,  MD  21010 

Commander/Di  rector 
Chemical  Systems  Laboratory 
US  Arny  Armament  Research 
4  Development  Command 
ATTN:  DRDAR-CLB-PS  (Mr.  Joseph  Vervier) 
Aberdeen  Proving  Ground,  MD  21010 

Commander/Di  rector 
Chemical  Systems  Laboratory 
US  Army  Armament  Research 
&  Development  Command 
ATTN:  DRDAR-CLY-A  (Mr.  Ronnald  Pennsyle) 
Aberdeen  Proving  Ground,  MD  21010 

Commander 

US  Army  Ballistic  Research  Laboratory/ 
ARRADCOM 

ATTN:  DRDAR-TSB-S  (STINFO) 

Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Army  Electronics  Research 
&  Development  Command 
ATTN:  DRDEL-CCM  (W.  H.  Pepper) 

Adel  phi,  MD  20783 

Commander 

US  Arny  Electronics  Research 
&  Development  Command 
ATTN:  DRDEL-CG/DRDEL-DC/DRDEL-CS 
2800  Powder  Mill  Road 
Adel  phi,  MD  20783 


Commander 

US  Army  Electronics  Research 
&  Development  Command 
ATTN:  DRDEL-CT 
2800  Powder  Min  Road 
Adel  phi,  MD  20783 

Commander 

US  Army  Electronics  Research 
&  Development  Command 
ATTN:  DRDEL-PAO  (Mr.  Steven  Klmmel ) 

2800  Powder  Mill  Road 
Adel  phi,  MD  20783 

Project  Manager 
Smoke/Obscurants 
ATTN:  DRDPM-SMK 

(Dr.  Anthony  Van  de  Wal ,  Jr.) 
Aberdeen  Proving  Ground,  MD  21005 

Project  Manager 
Smoke/Obscurants 

ATTN:  DRDPM-SMK -T  (Mr.  Sidney  Gerard) 
Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Army  Test  A  Evaluation  Command 
ATTN:  DRSTE-AD-M  (Mr.  Warren  M.  Baity) 
Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Army  Test  &  Evaluation  Command 
ATTN:  DRSTE-AD-M  (Dr.  Norman  E.  Pentz) 
Aberdeen  Proving  Ground,  MD  21005 

Director 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-AAM  (Mr.  William  Smith) 
Aberdeen  Proving  Ground,  MD  21005 

Di rector 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-CS  (Mr.  Philip  H.  Beavers) 
Aberdeen  Proving  Ground,  MD  21005 

Director 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-GB  (Wilbur  L.  Warfield) 
Aberdeen  Proving  Ground,  MD  21005 


Director 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-GP  (Mr.  Fred  Campbel 1 ) 

Aberdeen  Proving  Ground,  MD  21005 

Di rector 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-GS 

(Mr.  Michael  Starks/Mr.  Julian  Chernlck) 
Aberdeen  Proving  Ground,  MD  21005 

Di rector 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-J  (Mr  James  F.  O'Bryon) 
Aberdeen  Proving  Ground,  MD  21005 

Di rector 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-LM  (Mr.  Robert  M.  Marchetti ) 
Aberdeen  Proving  Ground,  MD  21005 

Commander 

Harry  Diamond  Laboratories 
ATTN:  Dr.  William  W.  Carter 
2800  Powder  Mill  Road 
Adel  phi,  MD  20783 

Commander 

Harry  Diamond  Laboratories 

ATTN:  DELHD-R-CM  (Mr.  Robert  McCoskey) 

2800  Powder  Mill  Road 
Adel  phi,  MD  20783 

Commander 

Harry  Diamond  Laboratories 
ATTN:  DELHD-R-CM-NM  (Dr.  Robert  Humphrey) 
2800  Powder  Mill  Road 
Adel  phi,  MD  20783 

Commander 

Harry  Diamond  Laboratories 
ATTN:  DELHD-R-CM-NM  (Dr.  Z.  G.  Sztankay) 
2800  Powder  Mill  Road 
Adel  phi,  MD  20783 

Commander 

Harry  Diamond  Laboratories 
ATTN:  DELHD-R-CM-NM  (Dr.  Joseph  Nemarlch) 
2800  Powder  Mill  Road 
Adel  phi,  MD  20783 
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Commander 

Air  Force  Systems  Command 

ATTN:  WER  (Mr.  Richard  F.  Plcanso) 

Andrews  AFB,  MO  20334 

Martin  Marietta  Laboratories 
ATTN:  Jar  Mo  Chen 
1450  South  Rolling  Road 
Baltimore,  MO  21227 

Commander 

US  Army  Concepts  Analysis  Agency 
ATTN:  CSCA-SMC  (Mr.  Hal  E.  Hock) 

8120  Woodmont  Avenue 
Bethesda,  MD  20014 

Director 

National  Security  Agency 
ATTN:  R52/Dr.  Douglas  Woods 
Fort  George  G.  Meade,  MD  20755 

Chief 

Intelligence  Materiel  Development 
A  Support  Office 

US  Army  Electronic  Warfare  Laboratory 
ATTN:  DELEW- I  (LTC  Kenneth  E.  Thomas) 
Fort  George  G.  Meade,  MD  20755 

The  John  Hopkins  University 
Applied  Physics  Laboratory 
ATTN:  Dr.  Michael  J.  Lun 
John  Hopkins  Road 
Laurel 1,  MD  20810 

Dr.  Stephen  T.  Hanley 
1720  Rhodesia  Avenue 
Oxon  Hill ,  MD  20022  * 

Science  Applications  Inc. 

ATTN:  Mr.  G.  D.  Currie 
15  Research  Drive 
Ann  Arbor,  MI  48103 

Science  Applications  Inc. 

ATTN:  Dr.  Robert  E.  Turner 
15  Research  Drive 
Ann  Arbor,  MI  48103 


Commander 

US  Army  Tank -Automotive  Research 
&  Development  Command 
ATTN:  DRDTA-ZSC  (Mr.  Harry  Young) 

Warren,  MI  48090 

Commander 

US  Army  Tank  Automotive  Research 
A  Development  Command 
ATTN:  DRDTA-ZSC  (Mr.  Wallace  Mick,  Jr.) 
Warren,  MI  48090 

Dr.  A.  D.  Belmont 
Research  Division 
Control  Data  Corporation 
P.0.  Box  1249 
Minneapolis,  MN  55440 

Director 

US  Army  Engr  Waterways  Experiment  Station 
ATTN:  WESEN  (Mr.  James  Mason) 

P.0.  Box  631 
Vicksburg,  MS  39180 

Commander 

US  Army  Research  Office 
ATTN:  DRXRO-GS  (Dr.  Leo  Alpert) 

P.0.  Box  12211 

Research  Triangle  Park,  NC  27709 
Commander 

US  Army  Research  Office 
ATTN:  DRXRO-PP  (Brenda  Mann) 

P.0.  Box  12211 

Research  Triangle  Park,  NC  27709 
Commander 

US  Army  Cold  Regions  Research 
A  Engineering  Laboratory 
ATTN:  CRREL-RD  (Dr.  K.  F.  Sterrett) 
Hanover,  NH  03755 

Commander/Di  rector 
US  Army  Cold  Regions  Research 
A  Engineering  Laboratory 
ATTN:  CRREL-RG  (Mr.  George  Altken) 
Hanover,  NH  03755 
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Commander 

US  Army  Cold  Regions  Research 
&  Engineering  Laboratory 
ATTN:  CRREL-RG  (Mr.  Roger  H.  Berger) 
Hanover,  NH  03755 

Commander 

US  Army  Armament  Research 
&  Development  Command 
ATTN:  DRDAR-AC  (Mr.  James  Greenfield) 
Dover,  NJ  07801 

Commander 

US  Army  Armament  Research 
&  Development  Command 
ATTN:  DRDAR-TSS  (Bldg  #59) 

Dover,  NJ  07801 

Commander 

US  Army  Armament  Research 
&  Development  Command 
ATTN:  ORCPM-CAWS-EI  (Mr.  Peteris  Jansons) 
Dover,  NJ  07801 

Command er 

US  Army  Armament  Research 
A  Development  Command 
ATTN:  DRCPM-CAWS-EI  (Mr.  G.  H.  Waldron) 
Dover,  NJ  07801 

Deputy  Joint  Project  Manager 
for  Navy/USMC  SAL  GP 

ATTN:  DRCPM-CAWS-NV  (CPT  Joseph  Micell) 
Dover,  NJ  07801 

Commander/Di  rector 
US  Army  Combat  Surveillance  A  Target 
Acquisition  Laboratory 
ATTN:  DELCS-I  (Mr.  David  Longlnotti) 

Fort  Monmouth,  NJ  07703 

Commander/Di  rector 
US  Army  Combat  Surveillance  A  Target 
Acquisition  Laboratory 
ATTN:  OELCS-PE  (Mr.  Ben  A.  01  Campli) 

Fort  Monmouth,  NJ  07703 


Commander/Di  rector 

US  Army  Combat  Surveillance  A  Target 
Acquisition  Laboratory 
ATTN:  DELCS-R-S  (Mr.  Donald  L.  Foiani) 
Fort  Monmouth,  MJ  07703 

D1 rector 

US  Army  Electronics  Technology  A 
Devices  Laboratory 
ATTN:  DELET-DD  (S.  Danko) 

Fort  Monmouth,  NJ  07703 

Project  Manager 
FIREFINDER/REMBASS 

ATTN:  DRCPM-FFR-TM  (Mr.  John  M.  Bialo) 
Fort  Monmouth,  NJ  07703 

Commander 

US  Army  Electronics  Research 
A  Development  Command 
ATTN:  DRDEL-SA  (Dr.  Walter  S.  McAfee) 
Fort  Monmouth,  NJ  07703 

OLA,  2WS  (MAC) 

Holloman  AFB,  NM  88330 

Commander 

Air  Force  Weapons  Laboratory 
ATTN:  AFWL/WE  (MAJ  John  R.  El  rick) 
Kirtland,  AFB,  NM  87117 

Di rector 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-SL 

White  Sands  Missile  Range,  NM  88002 
Director 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-SL  (Dolores  Angulano) 

White  Sands  Missile  Range,  NM  88002 

Director 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-TDB  (Mr.  Louie  Dominguez) 
White  Sands  Missile  Range,  NM  88002 

Di rector 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-TDB  (Mr.  William  J.  Leach) 
White  Sands  Missile  Range,  NM  88002 
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Director 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-TGP  (Mr.  Roger  F.  Willis) 

White  Sands  Missile  Range,  NM  88002 

Di rector 

Office  of  Missile  Electronic  Warfare 
ATTN:  DELEW-M-STO  (Dr.  Steven  Kovel) 

White  Sands  Missile  Range,  NM  88002 

Office  of  the  Test  Director 
Joint  Services  EO  GW  CM  Test  Program 
ATTN:  DRXDE-TD  (Mr.  Weldon  Findley) 

White  Sands  Missile  Range,  NM  88002 

Commander 

US  Army  White  Sands  Missile  Range 
ATTN:  STEWS-PT-AL  (Laurel  B.  Saunders) 
White  Sands  Missile  Range,  NM  88002 

Commander 

US  Army  R&D  Coordinator 
US  Embassy  -  Bonn 
Box  165 

APO  New  York  09080 

Grumman  Aerospace  Corporation 
Research  Department  -  MS  A08-35 
ATTN:  John  E.  A.  Selby 
Bethpage,  NY  11714 

Rome  Air  Development  Center 
ATTN:  Documents  Library 
TSLD  (Bette  Smith) 

Griffiss  AFB,  NY  13441 

Dr.  Roberto  Vagllo-Laurin 
Faculty  of  Arts  and  Science 
Dept,  of  Applied  Science 
26-36  Stuyvesant  Street 
New  York,  NY  10003 

Air  Force  Wright  Aeronautical  Laboratories/ 
Avionics  Laboratory 

ATTN:  AFWAL/AARI-3  (Mr.  Harold  Geltmacher) 
Wright-Patterson  AFB,  OH  45433 


Air  Force  Wright  Aeronautical  Laboratories/ 
Avionics  Laboratory 

ATTN:  AFWAL/AARI-3  (CPT  William  C.  Smith) 
Wright-Patterson  AFB,  OH  45433 

Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSF-CF-R  (CPT  James  M.  Watson) 

Fort  Sill,  OK  73503 

Commandant 

US  Army>F1eld  Artillery  School 
ATTN:  ATSF-CD-MS 
Fort  Sill,  OK  73503 


Commandant 

US  Army  Field  Artillery  School 
ATTN:  NOAA  Liaison  Officer 

(CDR  Jeffrey  G.  Carlen) 

Fort  Sill,  OK  73503 

Commandant 

US  Army  Field  Artillery  School 
Morris  Swett  Library 
ATTN:  Reference  Librarian 
Fort  Sill,  OK  73503 

Commander 

Naval  Air  Development  Center 
ATTN:  Code  301  (Mr.  George  F.  Eck) 
Warminster,  PA  18974 

The  University  of  Texas  at  El  Paso 
Electrical  Engineering  Department 
ATTN:  Dr.  Joseph  H.  Plerluissi 
El  Paso,  TX  79968 

Commandant 

US  Amy  Air  Defense  School 

ATTN:  ATSA-CD-SC-A  (CPT  Charles  T.  Thorn) 

Fort  Bliss,  TX  79916 


Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSF-CF-R 
Fort  Sill,  OK  73503 
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Commander 

HQ,  TRADOC  Cobmlned  Arms  Test  Activity 
ATTN:  ATCAT-OP-Q  (CPT  Henry  C.  Cobb,  Jr.) 
Fort  Hood,  TX  76544 

Commander 

HQ,  TRADOC  Combined  Arms  Test  Activity 
ATTN:  ATCAT-SCI  (Dr.  Oarrell  W.  Collier) 
Fort  Hood,  TX  76544 

Commander 

US  Army  Dugway  Proving  Ground 
ATTN:  STEDP-MT-DA-L 
Dugway ,  UT  84022 

Commander 

US  Army  Dugway  Proving  Ground 

ATTN:  STEDP-MT-DA-M  (Mr.  Paul  E.  Carlson) 

Dugway,  UT  84022 

Coipander 

US  Army  Dugway  Proving  Ground 

ATTN:  STEDP-MT-DA-T  (Mr.  John  Trethewey) 

Dugway,  UT  84022 

Commander 

US  Army  Dugway  Proving  Ground 

ATTN:  STEDP-MT-DA-T  (Mr.  William  Peterson) 

Dugway ,  UT  84022 

Defense  Documentation  Center 
ATTN:  DDC-TCA 
Cameron  Station  Bldg  5 
Alexandria,  VA  22314 
12 

Ballistic  Missile  Defense  Program  Office 
ATTN:  DACS-BMT  (Colonel  Harry  F.  Ennis) 
5001  Elsenhower  Avenue 
Alexandria,  VA  22333 

Defense  Technical  Information  Center 
ATTN:  DDA-2  (Mr.  James  E.  Shafer) 

Cameron  Station,  Bldg  5 
Alexandria,  VA  22314 


Commander 

US  Amy  Materiel  Development 
A  Readiness  Command 

ATTN:  DRCBSI-EE  (Mr.  Albert  Glambalvo) 
5001  Elsenhower  Avenue 
Alexandria,  VA  22333 

Commander 

US  Army  Materiel  Development 
A  Readiness  Command 
ATTN:  DRCLDC  (Mr.  James  Bender) 

5001  Elsenhower  Avenue 
Alexandria,  VA  22333 

Defense  Advanced  Rsch  Projects  Agency 
ATTN:  Steve  Zakanyez 
1400  Wilson  Blvd 
Arlington,  VA  22209 

Defense  Advanced  Rsch  Projects  Agency 
ATTN:  Dr.  James  Tegnelia 
1400  Wilson  Blvd 
Arlington,  VA  22209 

Institute  for  Defense  Analyses 
ATTN:  Mr.  Luclen  M.  Blberman 
400  Army -Navy  Drive 
Arlington,  VA  22202 

Institute  for  Defense  Analyses 
ATTN:  Dr.  Ernest  Bauer 
400  Army-Navy  Drive 
Arlington,  VA  22202 

Institute  of  Defense  Analyses 
ATTN:  Dr.  Hans  G.  Wolfhard 
400  Army -Navy  Drive 
Arllngon,  VA  22202 

System  Planning  Corporation 
ATTN:  Mr.  Daniel  Friedman 
1500  Wilson  Boulevard 
Arlington,  VA  22209 

System  Planning  Corporation 
ATTN:  COL  Hank  Shelton 
1500  Wilson  Boulevard 
Arlington,  VA  22209 
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US  Army  Intelligence  4  Security  Command 
ATTN:  Edwin  Speakman,  Scientific  Advisor 
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